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ABSTRACT 

Consider binary system with a millisecond pulsar ejecting relativistic particles and an 
optical star emitting soft photons with the energy u> ~ 1 — 10 eV. These low-energy 
photons are scattered by the relativistic electrons and positrons of the pulsar's wind. 
The scattered photons forms a wide spectrum from hard X-ray band up to gamma 
band e ~ 1 — lOOOGeV .When the pulsar wind is isotropic the luminosity of gamma 
radiation L 7 = Ly(ip) depends heavily on the angle tp between the directions to the 
optical star and to the observer from the pulsar. During the orbital motion this angle 
varies periodically giving rise to the periodical change of the observed intensity of the 
gamma radiation and its spectrum. We calculated the spectral shape of the scattered 
hard photons. Under the assumption that the energy losses of the relativistic particles 
are small we receive analytical formulas. We apply our results to the binary system 
PSR B1259-63 and show that if the wind from the Be star is accounted for then it is 
possible to reproduce the observed spectrum. 

Key words: binaries: general - pulsars: general - pulsars:individual:PSR B1259-63 - 
X-rays:stars 



1 INTRODUCTION 

The principal accessible store of energy in a pulsar is its rota- 
tional energy, which is liberated at a rate L p ~ I V 4n 2 /p 3 ~ 
, 10 32 -r 10 36 erg/s ( see, for example, a table of pulsars pa- 
rameters ( L p ,p, P) given in the book of Beskin et al.,1993 ), 
where I ~ 10 45 g cm -2 is moment of inertia, p is a period, P is 
a deceleration of a pulsar. The bulk of the pulsar energy L p 
is transferred to the pulsar wind which consists of electrons, 
positrons and probably heavy ions, L w — gL p , g < 1. The 
Lorentz factor of the relativistic particles in the wind may 
vary m range 7 ~ 10- 10 6 (e.g. Manchester&Taylor,1977). 
The intrinsic gamma ray luminosity of pulsed emission from 
the short periodic pulsars is of the order of L^/L p w 0.01 
(see Arons, 1991). Below we discuss the different (induced) 
mechanism of non pulsed gamma radiation generation in 
binary with a pulsar and an optical star. This mechanism 
could result in considerably higher ratio of L 7 /L p . The pre- 
liminary results of this work were published in the paper of 
Chernyakova&Illarionov, 1997. 

Consider the case of binary with a pulsar ejecting rel- 
ativistic particles and an optical star emitting soft photons 
in optic and UV band with the energy ui ~ 1 — 10 eV. These 
low-energy photons are scattered by the pulsar wind rela- 
tivistic electrons and positrons. The energy of the photon af- 
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ter the inverse Compton scattering is very high - e max ~ <^7 2 
in the Thomson limit (luj <C mc 2 ) and e max ~ mc 2r y in the 
opposite case ( m is a mass of an electron, c is a light ve- 
locity). The scattered photons form a wide spectrum from 
hard X-ray band up to gamma band e ~ 1 — lOOOGeV. The 
relativistic particle scatters the photon preferably along the 
direction of the particle velocity. As a result while soft pho- 
tons are directed radially from the optical star, the scattered 
hard photons move radially from the pulsar. Here and bel- 
low we assume that the pulsar wind particles are radially 
directed from the pulsar inside the effective scattering vol- 
ume. 

In the case of the presence of an obstacle for the pulsar 
wind and mainly when the matter flow from the optical star 
is rather intensive the radial flow of the relativistic wind is 
destroyed. Then the system of shock waves resulting from 
the collision of two winds appears between the pulsar and 
the optical star and the trajectories of the electrons and 
positrons beyond the shock change. In the work of Tavani 
and Brookshaw (1991) the case of weak (in comparison with 
the pulsar wind) matter outflow is discussed. The hydrody- 
namics of the collision between the relativistic and the non- 
relativistic winds is closely analogous to the hydrodynamics 
of the two nonrelativistic winds collision, which was inten- 
sively discussed by different authors applied to the WR+OB 
binaries since the work of Prilutskii&Usov (1976). 

The total luminosity of scattered hard photons L-, is 
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equal to the total particle energy losses Li oss through the 
inverse Compton scattering in the course of the particle 
motion from the pulsar to the infinity. To estimate L 7 
we note that the rate of the energy losses of a relativis- 
tic particle moving in a radiation field with energy density 
of soft photons w so f t ~ L»/47ra 2 c is about mc^dj/dt ~ 
—WsoftCTTC-y 2 in the Thomson limit. Here L, is a luminos- 
ity of the optical star, a is a binary stars separation and 
ot = |(e 2 /mc 2 ) 2 = 6.65 x l(T 25 cm 2 is the Thomson 
cross-section. Hence the decrease of the Lorentz factor of a 
particle is A7 ~ 7 [l — (1 + 7/7*) -1 ] , 7* = 4iramc 3 /gtL-* . 
So the rate of all particles total energy loss is 

Lj ~ Lloss ~ ~T~ ~ V ~ 1 + 7/% ) ' 

Thus in the case 7^7, the gamma-ray luminosity from 
binary is proportional to the luminosity of the optical star 
L* and to the luminosity of the wind L w . L<y = KL W , where 
the transformation parameter K = 7/7*. In the close binary 
system with an optical star with high luminosity 7 7, . In 
this case practically all energy of the wind transfers to the 
energy of scattered photons L 

-7 max ~ L w . 

In each unit of the volume in the region where 
the gamma radiation is generated the source function of 
gamma radiation is highly anisotropic. When the pulsar 
wind is isotropic the gamma luminosity L 7 (V>) ( L-y = 
/ 2jrL 1 (ip)d cos ip) has an azimuthal symmetry around the 
line connected binary companions and depends heavily on 
the angle ip between the directions to the optical star and 
to the observer from the pulsar, see Figure 1. In our case 
of free radial relativistic wind the binary system emits the 
maximum energy of gamma radiation in the direction of the 
star (ip = 0) and the minimum energy of radiation in the 
opposite direction (ip = n). The spectrum of the radiation 
and the maximal radiated energy also depend on ip. Dur- 
ing orbital motion ip varies periodically giving rise to the 
periodical change of the intensity of the gamma radiation 
coming from the binary system to the observer. 

In section 2 we calculate the spectral shape L 1 (e,ip) of 
the scattered hard photons in the case of arbitrary value of 
the parameter -^z, going beyond the Thomson limit. Under 
the assumption ]f < 1 we receive analytical formula for the 
L 7 (e,V0- 

In section 3 we apply our results to the binary sys- 
tem PSR B1259-63 and find that under the assumption 
of the power law relativistic particles spectrum, d ° ± = 
0AL w /[mc 2 ( 7ml ° n 4 - 7m ° x 4 )] 7 - 2 ' 4 in the range 10 < 7 < 
500, our model describes the observe photon spectrum 
rather good but the intensity is less then the observed one by 
a factor about 30. This discrepancy is due to the presence of 
the mass outflow from the Be star which disturbed the free 
flow of the pulsar wind. The centrally located shock appears 
between the pulsar and the star due to the interaction be- 
tween the two winds. The big differences between the values 
of the velocities of the particles from the different sides of 
the tangential discontinuity will lead to the growth of the in- 
stabilities and the two winds will be macroscopically mixed 
between the shocks. Then the heavy non relativistic wind 
slows down the volumes filled by the relativistic electrons 
and positrons and they acquire essentially non relativistic 
hydrodynamic drift velocity va along the shock while the 



energy of electrons and positrons does not changes signifi- 
cantly. With the decrease of the hydrodynamic velocity of 
the relativistic plasma the time which it spends near the 
optical star increases in c/vd times. The effective transfor- 
mation parameter K e fj ~ — K thus can be large enough 
to overcome the discrepancy between the simple theory and 
observations. 



2 GENERATION OF GAMMA RADIATION 

Lets us consider an interaction between the relativistic pul- 
sar wind and the soft radiation from the companion and 
find the spectral and angular dependence of the outgoing 
hard radiation Lj(e,ip). We assume that the pulsar wind 
and the soft emission from the companion are isotropic and 
that there is no mass outflow from the optical star. We treat 
both the pulsar and the optical star as a point sources. In 
this case the trajectories of the particles and soft non scat- 
tered photons are directed radially from the pulsar and from 
the optical star correspondingly. 

Let I denote the location of the small element of volume 
dV at a distance r to the pulsar and at an angle 82 to the 
line of sight P(3 (see Figure 1). In this element of volume 
optical photons with the energy u> scatters by the relativis- 
tic particles (electrons or positrons). The angle between the 
directions of the photon and the particle movements before 
the interaction is denoted by 6\. The photon scattering an- 
gle is designated by 0. Lets I signifies the distance from this 
volume to the companion S. 

We are interested in the case of a distant observer and 
thus a vector along the direction to the observer from any 
point near the system may be considered as parallel to such 
a vector from the pulsar. Only the photons scattered in the 
direction of the observer within the small solid angle fi = 
S cosx/D 2 will reach the observer. Here S is the area of 
the observer's surface, D is a distance from the pulsar to 
the observer and \ is an angle between the direction to the 
observer from the pulsar and from the volume. As \ is of 
the order of a/D <C 1, we may neglect the variance of Q, 
from point to point. As we have already mention relativistic 
particle scatters the photons preferably along the direction 
of the particle velocity and thus 82 <C 1. From this fact it 
follows that angles 6 and 61 depe nd only on r and with an 
accuracy of the order 62 ~ \J ujJs are equal. 

2.1 Anisotropy of the energy radiated from the 
binary system in the Thomson limit. 

Lets calculate the ^-dependence of the luminosity of the 
hard radiation emitted from the binary system L y (ip) in the 
Thomson limit ( coy / mc 2 <C 1 ) . The luminosity of scattered 
hard photons L 7 (ip) is equal to the particle energy losses in 
the course of the particle motion from the pulsar to the in- 
finity. So at first lets calculate the total energy losses Li oss t 
of the electron moving along the radial trajectory at an an- 
gle ip to the line PS in the radiation field of the optical star. 
The rate of the relativistic ( 7 >• 1 ) electron's energy losses 
is given by 

1 d'y _ d'y L^ary 2 (1 — cos 81) 2 
c dt dr ATvmc 3 P 
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where 

I = W \r — a cos if)) 2 + a 2 sin 2 ip, 

cos 61 = (r — a cos ip)/ 1. 

Upon integrating ([j]) over radius from r — to r we 
obtain the dependence of the Lorentz factor of the electron 
on r : 



7o [1 + K$(i>)+ 
K / 2 



where 



$(V) = 



sini/) V VT+ 



3(tt - -0) 
2 sin i/i 



= arccotx 

,2 2 



2(1+: 



,(2) 



cos V 1 



(3) 



is a beaming function which as it will follows from ( 
represents the anisotropy of the gamma radiation, 

K = 70/7* = o"t£'*7o /(47ramc 3 ) 



(4) 



represents the efficiency of the electron energy transforma- 
tion to the energy of gamma radiation, the parameter x is 
determined as 

r/a — cos?/> 
sin ip 

and 70 is the initial Lorentz factor of the electron. Letting 
r^oowe find from (Q) the total change of the Lorentz fac- 
tor A7O) = 70 - j(x = 00, tp) = 70 [1 - (1 + ^(vor 1 ] • 

As the bulk of the energy losses of the relativistic elec- 
tron transfers to the scattered photons moving in the direc- 
tion of the electron movement the angular dependence of the 
energy radiated from the binary system Ly(ip) agree with 
the angular dependence of the energy losses of the electrons. 
The total energy losses of the isotropic flow of all relativis- 
tic particles in the direction of the observer located at an 
angle ip to the line PS in a small solid angle il are propor- 
tional to the amount of particles moving in this direction 
N e± Q/Aiv = L m r2/47r7omc 2 in the case of monoenergetic 
pulsar wind, and to the total energy losses of one particle 
mc 2 Ay(ip) 



L 7 (V>)« 



LlosaT[m= hE$m n 



4tt 



47T70 

1 



l + K<f>(ip) 



Q. 



(5) 



This is the energy transferred to the gamma band. Figure 2 
shows the dependence <&(?/>). The beaming function <&(ip) is 
rapidly decreasing function. At t/> < 1 <&(ip) — ir/ip and at 
n-ip < 1 $(?/>) ~ (7r--i/)) 4 /120. Equation (|j) is valid for the 
arbitrary value of K but it has no use for ip less then R* /a 
and bigger then n — R,/a (_R* - is a radius of the optical 
star) as the optical star is not a point source. That's why 
we don't worry about the fact that $(?/>) — > 00 as ip — > 0. 
If ip is not too small then it follows from (|^) that under the 
condition A'< 1 the energy decrease is small and 



itl)70 
4-7T7* 



(6) 



Upon integrating (|6|) over ip from ip — to ip = tv we obtain 
the total energy transferred to the gamma band in the case 
if< 1 



(§•■-») 



7oi« 



(7) 



2.2 The spectral and the angular dependence of 
the outgoing radiation 

Now lets calculate the spectral and the angular dependence 
of the outgoing radiation i 7 (w, e, 7, ip) in the case of ar- 
bitrary value of the parameter ujj/mc 2 , going beyond the 
Thomson limit. The luminosity of the scattered quanta mov- 
ing at an angle ip to the line PS within a unit solid angle in 
a unit of time in the energy range from e to e + de may be 
written in the following form: 



L~f (w, e, 7, ip) de = / en e ± n^c (1 — /3 cos 0\) ondV, 



(8) 



where j3 = v/c, n e ± is the density of the particles, 
nuj = L* I (4:ircl 2 Lu) is the density of the optical photons, 
<jk (0, 0i, 02 u) - is the Klein-Nishina cross-section (Jauch 
&Rohrlich 1976) 



3<tt s , a ,_2 
(1-/3 cos 0i ) x 



167T7 2 LU 2 



(9) 



7 2 (1 - (3 cos 61) (1 - /3cos0 2 ) 



- 1 



+ 1 



+ 



( 1 — cos ( 



mc 2 7 I (1 -/3 cos 02 )(!-/? cos 6>i) 



The energy of the scattered photon is ( Berestetskii et 
al. 1971 ) 



P cos ( 



1 -/?cos6> 2 + -V (1 



From ( |Io| ) it follows that 



cos t>2 



P- 1 



-(1 - /3 cos 6»i ) H -5- (1 - cosf 

e mc 2 -y 



(10) 



.(11) 



We are interesting in case 7 2> 1 and in all formulas that fol- 
lowed we use the expansion P = 1 — . In the spherical co- 
ordinates dV = 2irr 2 drd cos 02 = 2nr 2 uj ( 1 — cos #1 ) dedr / e 2 . 
When e S> w, then 02 <C 1 and thus 

r/a — costp 



cos Pi 



y (r/a — cos?/;) 2 + sin 2 ip 
Now it is possible to rewrite (^|) as an integral over r 
L-y (w,e,7, = 

r max 

AT W .A7,± u f „ „ x2 dr 



(12) 



87ra 2 c £ 



(1 — cos #i) ax 



(r/a— cos?/;) + sin 2 ?/> ' 



(13) 



where AT W is the number of photons going from the compan- 
ion per second and N e ± is the number of particles leaving 
the pulsar per second, m„ = N u j 'AttcI 2 , n e ± = N e ± /Ancr 2 , 
r max is determined from the condition cos02 = 1. If Lorentz 
factor of a particle is constant then from (|ll[), (^) follows: 
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a ^cos ip + 



sin ip 



1 

. z 

2 V z 



(14) 



where 



M = 



4^7 



A = 1 



In the case if <C 1 in Thomson limit the Lorentz lac- 
tor ol the particles is practically constant in the effective 
scattering volume. In the case wy/mc 2 > 1 this state- 
ment is also correct. In this case il the mean free path 
of a particle is much bigger then the binary separation: 
Ifree = a k i uU L„ / inuoa 2 c ^> a (ok SuU is a full Klein-Nishina 
cross-section) then the influence of the field of soft photons 
on the spectrum and flux of the relativistic particles may 
be neglected. In the highly relativistic case (uoy/mc 2 ^> 1) 
<JK Jull ~ armc 2 /toy and we note that from K <C 1 it follows 

that lf ree 3> O. 

After a substitution of the expression ^ in (jl3|) and 
the following integration under the assumption that 7 is 
constant we obtain : 



L y [fJ, e, 7, ip) = r F (uo, e, 7, ip) . 



where the constant r 

LO < £ < £max 



o 2 
Z7 UJ 



7^ UJ 

— and in the range 

1 + COS Ip 



1 + 



1 + COS ip) 



the function F(u},s, 7, ip) is determined as follows: 



F = 



simp 



arctan z — tp) 



(15) 



(16) 



(17) 



8 



3^l 22 ~ 3z tan i +tan i 



and is equal to zero when e > e n 
the equation (|l7]) simplifies: 

F = 2^t. 
sin 



When iii « e « £ t 



In the Thomson limit e max = 27 2 oj(1 + cosi/>) and hence 



is similar with a radiation from the black body with tem- 
perature T, then the number of soft photons in the energy 
range from oj to uj + duo is : 



dNu, = 



15L* uj duo 



tt 4 T 4 e<"/ T - 1 ' 
Replacing N u in formula ( |l5| ) with d/V^and integrating over 



to we receive: 



F b b (e, 7, 



15r JV e ±L» eFbb (e,7,^) 



1 

y2 



y2 



uoduo 



( e «,/T_ 1 y 



(19) 



(20) 



1-A = 



equation (1171) becomes 



< 2^J-(l + cosi/)) <g 1. Under these conditions L 7 (uo,e,jo,ip) = rN^Nt 



Figure 3 shows the dependence of Fbb on e for different 
values of ip and 7. It can be seen that in the energy range 
where the difference of the cross-section from the Thomson 
one is significant the form of the spectrum changes consid- 
erably. Figure 3 also illustrates the dependence of the maxi- 
mum energy of the scattered photons £ max on 7 ( according 
to ( p^[ ) within the Thomson limit £ max = 2ujy 2 (l + cost/;) 
and in the opposite case ^J-(l + cosip) l> 1 the maximum 
energy £ max = 2mc 2 y ). The ip— dependence of the energy 
gone with the scattered photons in a unit of time in a unit 
solid angle can be seen in Figure 4. The energy is normalized 
to C = ±5f N e ±L^ 2 . 

During the orbital motion ip changes periodically. In 
Figure 5 the dependence of Fbb on e for different positions 
of the companions in the circular orbit and for different val- 
ues of 7 is shown. In this figure i is an inclination angle 
and ip is a true anomaly. The relation between ip and (p is: 
cost/; = sin i cos i^. Figure 6 shows the orbital angle ip de- 
pendence of the luminosity £-y(yj) going in the direction of 
the observer for the different values of inclination angles and 
Lorentz factors. 

Formulas ( |l4| ) - ( pp| ) were received under the assump- 
tions that Lorentz factor is constant. If we consider in the 
Thomson limit the dependence of 7on r according to formula 
(^) then from formula ( |l3| ) the luminosity of the scattered 
radiation is: 

e 



2 -F 7 (uj,e,y ,ip) ■ 



F T 



sin ip 
8 



[4 arctan z — 2ip- 



(18) 



where 

Fy (w,e,7o, 



(21) 



(22) 



3/./ : 



4 



o 3 n 2 , y>, 3^' 

2z — ,iz tan h tan — 



7o 2 



1+ [7- 2 (l-cosfl 2 ) _1 -l] 2 dr 



7 (r / a — cos ip) + sin ip 



4w7 2 



If we integrate equation ( |l5| ) with F = Ft over the energy 
of the scattered photons s we receive the ^-dependence of 
the total outgoing energy in according with the equation (|^). 
The relationship between functions <E> (ip) and Ft is: 



" ;F T - d£ 



uj-y^- 



UJf^ 



Under the real conditions the spectrum of the optical 
photons is not monochromatic. If the companion's radiation 



The dependences of r max and cos 82 on r are determined by 
the formulas @,(||). 

Figure 7 represents spectra of the scattered radiation 
obtained in the Thomson limit for different values of the 
parameter K = 70/7* with 70 = 10 4 . Spectra are calculated 
taking into account the dependence of 7 on r. Curve K = 
corresponds to the case of constant Lorentz factor, which is 
described by formula (^|). It can be seen from this figure 
that the bigger is the value of parameter K at a given angle, 
the softer is the spectrum of the hard radiation. This figure 
also illustrates that with the growth of parameter K the 
total energy of the outgoing radiation increases tending to 
L w as K — » 00. 



© 0000 RAS, MNRAS 000, 000-000 



Non-pulsed gamma radiation from binary system with a pulsar 5 



3 DISCUSSION 

We considered the scattering of the soft photons emitted by 
the optical star by the relativistic electrons and positrons 
from the pulsar wind. We showed that the intensity of the 
hard radiation came about from such a scattering is propor- 
tional to the luminosities of the pulsar and the optical star. 
But it is worth to mention that the luminosity of optical 
star itself depends on the luminosity of the pulsar as it ab- 
sorbs and reradiate part of the energy of the pulsar wind. 
If we denote by L st the intrinsic ( un-illuminated) luminos- 
ity of the star and by Li n d = /Z/ p i?*/4a 2 the energy of the 
pulsar reradiated by the optical star ( /is the fraction of 
the energy flux from the pulsar incident on the surface of 
the optical star that is converted into soft emission) then 
the total luminosity of the star is L to t = L s t + L in d- PSR 
1957-20 is an example of the system where the illumination 
of the optical star by a pulsar affects the luminosity of the 
optical star. A detailed study of the optical light curve of 
PSR 1957+20 done by Callanan et al. (1995) shows that 
the secondary is close to fill its Roche lobe with 7% - 20% 
of the incident flux converted to optical emission. Thus the 
described mechanism of the generation of the hard radiation 
is important not only in case of a highly luminous optical 
star in binary but also in case of a close binary with bright 
pulsar and an optical star with a big radius. 



3.1 Comparison with the observations. 

The only known binary system which contains radio pulsar 
and emits non pulsed X-ray radiation is PSR B1259-63 sys- 
tem. Since its discovery, the PSR B1259-63 system has been 
observed several times at X-ray energies. ROSAT observed 
the PSR B1259-63 system near to apastron in 1991-1992 
(Cominsky et al. 1994, Greiner et al. 1995). In January 1994 
the X-ray emission from the system PSR B 1259-63 during 
the periastron passage was observed by telescopes ASCA 
and OSSE (Grove et al. 1995, Kaspi et al. 1995). ASCA also 
observed the post-periastron emission from this system on 
February 28, 1994. 

The PSR B1259-63 system contains the radio pulsar 
with the spin period P = 47.76ms, rotating around the 
massive Be star SS 2883. According to Manchester et al., 
1995 the orbital eccentricity is e = 0.87, the projected 
semimajor axis is a mj sini = 3.9 x 10 13 cm and the lon- 
gitude of periastron is ui = 139°. Spin-down luminosity 
of the pulsar is L p ~ 9 x 10 35 erg/s. The mass function 
f(M p ) = ^"^y = 1-53M Q , the luminosity of the 

Be star is L* = 2.2 x 10 38 erg/s (Johnston et al. 1992, 
1994). The distance between companions at periastron is 
a, = a m j(l — e) ~ 10 13 cm. 

The review of the different discussed models of the ori- 
gin of the non pulsed X-ray spectrum from the system can be 
found in the paper of Tavani&Arons 1997. We tried to apply 
the results of our paper to this system. Assuming the power 
law distribution of the electrons and positrons in the pulsar 
wind 



Fav (e, i>) 



■Fw, (e, 7, •!/;) 7 " d 7- 



(24) 



dN e ± 

d-7 



Cy-°L W W, C = [(s - 2) /( 7 ^ n s - 7^) ] 
and rewriting formulas (|l9|), ( pcj ) taking into the account 
the dependence iV e ± (7)we receive 



The result of our model for 10 < 7 < 500, s = 2.4 and 
the measured spectrum are shown in Figure 8. Spectrum 
measured by OSSE (black dots) was taken from the paper of 
Grove et al. 1995. The straight lines shows the extrapolation 
of ASCA data (Kaspi et al. 1995). As it can be seen from the 
Figure 8 while the form of the spectrum is quite close to the 
observed one the intensity of the radiation resulted from our 
model is less then the observed one by a factor of 30. 

This discrepancy is due to the fact that we derive the 
formula ( p3[ ) under the assumption of the undisturbed free 
flow of the pulsar wind which is not correct in the system 
PSR B 1259-63 where the strong mass outflow from the Be 
star presents. Writing the mass rate of the polar wind of the 

Be star as M~ (10 -8 ' Mq / 'yr) M-sand the wind velocity as 
v = 10 8 «8cm/s we find the ratio of the impulses of the two 

winds Lp/cv M= 0.45/ns Ms ■ Thus for the Be star polar 
wind typical parameters the centrally located shock between 
the pulsar and the star due to the interaction between the 
two winds seems to appear and thus the free flow of the 
pulsar wind will be disturbed. The big differences between 
the values of the velocities of the particles from the different 
sides of the tangential discontinuity will lead to the growth 
of the instabilities and the two winds will be macroscopi- 
cally mixed between the shocks. The results of numerical 
calculations of Igumenshchev (1997) verified such a picture. 
Then the heavy non relativistic wind slows down the vol- 
umes filled by the relativistic electrons and positrons and 
they acquire essentially non relativistic hydrodynamic drift 
velocity along the shock while the energy of electrons and 
positrons does not changes significantly. With the decrease 
of the hydrodynamic velocity of the relativistic plasma the 
time which it spends near the optical star increases in c/vd 
times. The idea of small drift velocity as applied to the sys- 
tem LSI 61°303 was discussed by Maraschi&Treves (1981) 



-^K thus 



15r CL W L, eFav (e, ip) 
7T 4 mc 2 T 5 ' 



(23) 



The effective transformation parameter K e ff 
can be large enough to overcome the discrepancy between 
the simple theory and observations. 

Under the assumption of the power law energy distri- 
bution of the relativistic electrons and positrons after the 
shock d c± = A"f~ 3 and their isotropic velocity distribu- 
tion we can estimate the intensity of the X-ray radiation 
from the system. Such a distribution can be either the re- 
sult of the particle acceleration on the shock or in the case 
of the intrinsic power law distribution of the electrons and 
positrons in the pulsar wind and thin collisionless shock, 
passing which the relativistic particles do not change their 
energy distribution. 

Lets consider an element of volume dV filled with the 
relativistic plasma locating in the shock on a distance R 
from the Be star. In the case of the relativistic electrons 
and positrons isotropic distribution the number of particles 
produces the photons with energy e moving in the direc- 
tion of the observer in a unit solid angle is ^ ^dr/ = 

0.5 — gfr cicos #2^7 = \ dN d c± f? (1 — cos 8) ded'y (see section 
2.2). Then according to (j^, |], ]l^) the number of photons 
coming to the observer in the unit of time per unit square 
per MeV is 
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dN = — — x 

32tt luD 2 



-y/e/2^(l-cos 9) 



W7 2 (l — cos ( 



djdV(25) 



where 9 is a photon scattering angle. For the simplicity we 
take the photon density along the shock constant and equal 
to the one at a distance R = 10 13 cm. The factor AVc&n be 
estimated from the energy conservation law by equating the 
energy enters the shock per second L p £l/4nand the energy 
leaving the shock per second iVmc^^^/afs — 2), where 
f2 is the solid angle under which the shock wave is seen from 
the pulsar. Thus we have AV ~ L p Qa(s — 2)/47rmc 2 Vd7^~^. 
Then integrating ( ]25| ) over the volume of the shock we have 
for the N 

■■- L*L p Sl (s- 2) (2hO<-W a _- (1+ . )/a 

(l + «) ' ( ' 



N ~ 2r 



D 2 R 2 mcvd 7r 7 2 



For 7„ 



10, s = 2.4, Z? = 2kpc, v d = 10 8 cm/s we have 



AT ~ 5 x 10" 3 ( 



lOOfceV" 



7 ph/s/cm 2 /MeV 



(27) 



while the observable value of the radiation at a 100 keV is 
2.8 x 10~ 3 ph/s/cm 2 /MeV. Thus if the small drift velosity 
is taken into account it is possible to explain the observed 
spectral shape and intensity of the X-ray radiation. It is 
also worth to mention that for the particles with big Lorentz 
factor 7 > 7*1^/0 the assumption of constant Lorentz factor 
will be not valid due to the inverse Compton losses and the 
break in the photon spectrum will appear. The break in the 
photon spectrum at e — uj^Iv^/c? will appear. The index of 
the photon spectrum after the break will be bi gge r then the 
original one at 1/2. It can be also seen from (E5|) that the 
focus effect will take place - the intensity of the radiation 
is proportional to the (1 — cos S)' 1+s - )/ ' 2 and thus the major 
part of the radiation will be emitted toward the direction of 
the optical star. 
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Figure Captions 

1) The geometry under consideration. P is a pulsar, S 
is an optical star, / is a point of interaction of the electrons 
and photons, O is an observer. Note that these points are 
not in one plane in general case. 

2) The dependence of the beaming function $ on an an- 
gle tp. 

3) The dependence of Fbbon e for different values of ip 
for 7 = 10 4 and 7 = 10 6 . 

4) The -^-dependence of the energy gone with the scat- 
tered photons in a unit of time in a unit solid angle. The 
energy is normlized to C = ^ N e ± L r j 2 . 

5) The dependence of Ftton e for different positions of 
the companions in the circular orbit for 7 = 10 4 and 7 = 10 6 . 

6) The orbital angle tp dependence of the total energy 
gone with the scattered quanta in the unit of time in the 
unit solid angle for the different values of inclination angles 
a) 7 = 10 4 , b) 7 = 10 6 . 

7) Spectra of the scattered radiation obtained for differ- 
ent values of the parameter K = 70 /7* . Spectra are calcu- 
lated taking into account the dependence of 7 on r. 

8) Wc apply our model to the system PSR B1259-63. 
Theoretical result (solid line), OSSE spectrum of emission 
from 1994 January 3-23(black dots) and schematic extrapo- 
lations to the power-law spectra of ASCA observations from 
1993 December 28(dashed line), 1994 January 10 (dotted 
line), and 1994 January 26 (dashed-dotted line) are shown. 
ASCA extrapolations obtained from the analysis of Kaspi et 
al. 1995, OSSE results are taken from the paper of Grove et 
al. 1995. 
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